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Abstract. — A novel atomic beam splitter, using reflection of atoms off an evanescent light 
wave, is investigated theoretically. The intensity or frequency of the light is modulated in order 
to create sidebands on the reflected de Broglie wave. The weights and phases of the various side- 
vands are calculated using three different approaches: the Born approximation, a semiclassical 
path integral approach, and a numerical solution of the time-dependent Schrdinger equation. 
We show how this modulated mirror could be used to build practical atomic interferometers. 



1. Introduction. 

A number of experimental techniques have been developed to enable the interference of atoms 
to be observed, as the present special issue illustrates. The essential requirement for producing 
quantum interference is that a system can pass between two points in its configuration space 
via more than one path — that is, the quantum amplitude for a passage via either path is 
non-negligible for a given evolution of the system. The two paths can be visualized as forming 
a closed loop in configuration-space. If the relative phase between the two quantum amplitudes 
can vary, then the system can finish in one of two (or more) final states, with probabilities 
depending on this phase. 

The well-known features just mentioned are illustrated by the three main elements of a 
typical particle interferometer: a first beam splitter, one or more mirrors to bring the two 
paths back together in position, and a second beam splitter to close the loop. Diffraction or 
refraction can also be used to play the role of a "mirror" in bending an otherwise straight path 
for the interfering particle. 

(*) Unite de recherche associee au CNRS. 
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In this paper we describe a new beam splitter for neutral atoms, and also show how a practical 

interferometer can be made, having a number of promising properties such as simplicity, and 
fairly high (about 8%) transmission efficiency into the useful output states. The beam splitter 
is a vibrating mirror — that is a surface which reflects atom incident upon it, and which 
moves rapidly to and fro along the normal direction. Such a beam splitter can form the basis 
of a number of interferometer designs. An especially interesting possibility is a very simple 
interferometer having only a single "optical element" — a horizontal mirror is used, and atoms 
bounce repeatedly on it. Gravity plays the essential role of bringing the atomic trajectories 
back to the mirror surface, and the same vibrating mirror is used to separate and recombine 
the interferometer arms. 

The ability to create a variety of motions of the mirror surface enables one to manipulate 
the reflected de Broglie waves in a general manner — both delicate adjustments and large 
shifts of the atomic momentum can be produced. In the general case, one notes that at normal 
incidence, the path length for a wave to travel along the ^-axis from a position z, be reflected 
at Zm , and then return to z, is 2{z — Zm)- By varying Zm in time, the variation in path length 
can be understood as forming an "optical element" such as a prism, lens, or phase grating. 
Here the "lens" (or other clement) is extended in time and affects the motion along the mirror 
normal z, while a conventional lens is extended in position and affects the motion transverse to 
its axis. Narrow slits and amplitude gratings can be made by switching the mirror reflectivity 
between one and zero. Our proposal relics on the possibility of vibrating the mirror rapidly 
(typical vibration frequencies are in the MHz region). This would be difficult for a traditional 
mirror made of matter, but is easy to achieve for a mirror formed by a light field, since laser 
beam intensities or frequencies can be modulated rapidly using acousto-optic modulators. 

In the following, we first briefly consider beam splitter in general, and the basic principle of 
the vibrating mirror (Sect. 2). An ideal mirror for atoms would consist of a sharp potential 
barrier [1]. In practice, one cannot flnd an ideal mirror, and an important type of mirror for 
atoms is a quasi- resonant evanescent light wave at the surface of a dielectric [2] . This produces 
a potential y exp(— 2kz) having an exponential dependence on atomic position z. We consider 
atoms bouncing on such a potential, and calculate the effect of a time-modulation of the 
amplitude V{t) — for the exponential function, such a modulation is equivalent to moving the 
potential along the z-axis by z^^ given by 



y(i)e-2«^ = F(0)e-2'^("-^»W) or 2KZm{t) =\n{V{t)/V{0)). (1.1) 

In sections 3 and 4 we consider two perturbative methods to calculate the probability for 
an atom to be scattered by the vibrating mirror from one energy eigenstate (or plane wave) 
to another. The first method uses the Born approximation in first order; this approximation 
is valid when the first order scattering probability is low. The second method calculates the 
phase accumulated by an atom undergoing refiection from a vibrating mirror, using an action 
integral along the classical path of an atom reflected by a stationary mirror. The range of 
validity now includes the physically interesting situation where the scattering probability from 
one eigenstate to another is of the order of 1. In section 5 we consider the case of an initial 
wave-packet rather than a single plane wave, and compare the approximate analytic results 
with those of a numerical solution of the Schrodinger equation. This enables us to confirm the 
validity of the various methods. In section 6 we then go on to consider the application of these 
ideas to make a realistic atom interferometer. 
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2. Atom beam splitters. 

Up to now many atomic beam splitters have been suggested, and several demonstrated. In 
the following list we mention the beam splitter used in existing intcrfcromcters. These are, 
to our knowledge, the longitudinal Stark method [3]; the simple Young's slits arrangement 
[4, 5]; micro-fabricated gratings [6]; the Raman pulse technique [7, 8]; the optical Ramsey 
interferometer [9, 10]; and the longitudinal Stern-Gerlach interferometer [11]. There has been 
rapid progress and some very promising results, in term of the minimum detectable phase 
shift in a given integration time, combined with a large effective area of the interferometer — 
both are important since for most experiments the phase change produced by the effects to be 
measured is proportional to the area of the loop in parameter-space. 

Methods to produce larger beam separations have been investigated, notably adiabatic pas- 
sage in a "dark" state [12, 13]; the magneto-optical beam splitter [14], and Bragg reflection 
at crystalline surfaces [15]. The Raman pulse method has also already been used to produce 
very large splittings by the use of many pulses. The vibrating mirror can produce a beam 
separation Ap of the order of several Kk for an incident momentum p = 100?i«;, as we will 
show. This allows a useful effective area for the interferometer without making it too sensitive 
to misalignments. (The parameter k in Eq. (1.1) is of the order of the wave vector for light 
in resonance with the atomic transition, so fin is approximately equal to the familiar "recoil" 
momentum.) 

The basic idea of the vibrating mirror is familiar from optics. To calculate the effect of a 
reflection from a mirror whose position varies sinusoidally, Zmit) = zq sin ivt, we assume the 
incident and refected waves can be written 

(t>inc{z,t) = expi{—kz — fit) (2.1) 

(f>iefi{z,t) ~ expi{kz — Clt — usinivt + it) (2.2) 

The reflected wave here is an approximate solution of the wave equation, valid when uoj <C CI. 
We look for a solution having a node on the mirror surface: 

MZm{t),t)=0 (2.3) 

This implies that 

u = 2kzo (2.4) 

The reflected wave has a carrier frequency CI plus a frequency modulation imposed by the 
mirror. It can be decomposed into its component frequencies as follows: 

oo 

ex.-pi{—Ctt — usiwoji) = ex])[—iCtt) J„(u) exp(— mwf). (2-5) 

n= — cx) 

The weight of a given sideband f2 ± nw is thus given by j J„(2fc2;o)P- 

Our vibrating mirror for matter waves works along the same general principles. An atom ar- 
riving with the momentum pi has an energy hCl = p?/2M. After reflection its final momentum 
Pf is given by 

2 2 

lL^ = Il^+nnuj (2.6) 
2M 2M ^ ' 

where n is a positive or negative integer. For hjo <C pf /2M, the momentum transfer Ap = pf — Pi 
is simply given by 

Ap ~ n = nq (2.7) 

Pi 
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where the elementary momentum spacing q is equal to 

q = 2.8 

Pi 

The efRciency of the transfer from pi to pf, for the case of a vibrating exponential potential, 
is the subject of the following sections. 

3. Perturbative calculation in the Born approximation. 

In this section we present a perturbative calculation of the momentum transfer for atoms re- 
flected off a modulated evanescent potential. Since the problem is invariant under a translation 
parallel to the mirror surface, we only consider the motion along the ^;-axis, defined normal to 

the surface. 

3.1. The coupling between unperturbed states. — The Hamiltonian is split into two 
parts: 

H = Ho + Vi 

where Hq is the un-modulated part: 

Ho = ^ + VoeM-'^i^z) (3.1) 

This Hamiltonian is responsible for the standard reflection of atoms off the evanescent wave 
[2]. Vb is the light shift of the ground state of the atom at the prism- vacuum interface located 
at z = 0. We consider that the atom adiabatically follows the energy state given by (3.1) 
without emitting any spontaneous photons. For simplicity, we further do not consider possible 
atom-surface interactions (van dcr Waals potential, e.g.), supposing that the atom remains far 
enough from the surface so that these are negligible compared to the evanescent wave potential 
in H. Vi is the modulated part of the atom-laser interaction [16]: 

Vi = eVo cxp(-2kz) sm{ujt) (3.2) 

We will treat the case that the modulation amplitude e is between and 1. 

We evaluate the efficiency of the momentum transfer along Oz by calculating perturbatively 
the coupling between cigcnstatcs of the unperturbed Hamiltonian Hq induced by the time 
dependent part Vi. The eigenstates of Hq, characterized by their asympotic momenta p > 0, 
are given by [17] (P = p/Uk is the scaled momentum): 

/op 

*p(^) = —smh{TTP)Kip[w{z)] (3.3) 
where Kip[w{z)] is the Bessel ii'-function of imaginary parameter iP, and 



w(zj = — exp[—KZ) (3.4) 

These wavefunctions are the eigenfunctions of the unperturbed Hamiltonian Hq, normalized 
in a box between z = and z = k~^. 
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According to Fermi's golden rule, the probability for an atom initially in '^'^{z) (correspond- 
ing to an incident momentum pi) to make a transition to ^{{z) (corresponding to a momentum 
pt) is given by: 

= l^|(vt'?|eyoexp(-2K^)|*P)|V(^f = i^i ± M (3.5) 
where p{Ef = Ei± Uto) is the density of states at the final energy: 

p{Ei = Ei±nw) = ^ = ^ (3.6) 
and $ corresponds to the flux of the incident atomic wave: 

The coupling term in (3.5) can be obtained analytically [18], yielding: 



^ sinh(|(Pi + Pf)) sinh(f (Pi - Pf)) 
We now discuss this result in the limits of large and small momenta Pi, Pf, respectively. 



^fi = «-h(.Pi) sinh(.Pf ) I o,^^ \ (3-8) 



3.2. Semiclassical limit. — In the limit Pi,Pf » 1, we may replace by exponentials the 
sinh function in (3.8) whose arguments are of the order of Pi,Pf . The transition probability is 
then approximated by (AP = Pf — Pi): 

The function 



appearing in (3.9) is equal to unity for x = and decreases exponentially for \x\ » 1. It 
describes the decrease of the efficiency of the momentum transfer if the potential undergoes 
too many oscillations during the reflection of the atom. Indeed the interaction time of the 
atom with the evanescent wave is given by 



^ (3.11) 

KPi 



and the momentum transfer can be written 



^..r.g (3.12) 

where Q = q/hn. The rapid decrease of the function P{x) means that only a few Hk of 
momentum can be transferred efficiently. 

In figure 1, we show Wa at pi = 100 Hk for the upper and the lower sideband, as a function 
of Q. Note that the results (3.8) and (3.9) predict an asymmetry in the transfer efficiency 
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eiemenlary momenlum spacing Q 



Fig. 1. — Heights of the sidebands for orders n = ±1 as a function of Q = lor. The incident 
momentum is pi = 100 Hk, and the modulation depth is e = 1. The values are normalized to an 
incident wave with unit amplitude. Dashed lines: the result of the Born approximation, Wa (Eq. 
(3.8)), multiplied by a factor pi/pt in order to normalize the amplitude of the incident wave. The 
upper (lower) line corresponds to n = +1( — 1). Solid line: result of the scmmiclassical approximation, 
|aip, (Eq. (4.20)). Circles: numerical solution (Sect. 5) for the normalized momentum distribution 
\'<p{p)/ip{pi)\^ with p^ — pf + 2Mhu> (i.e. sideband n — +1). For the sake of clarity, the points have 
been joined in the region Q < 4.2 (light dashed line). The parameters of the numerical calculation are 
given in the text (Sect. 5). Crosses: same as before, but now n = —1. 



between the upper and lower sideband Pf± = \/ ± 'iPiQ. This asymmetry is due to the fact 
that the monemtum transfer AP± = Pf± — Pj is not the same for the two sidebands (^): 

|AP+|-|AP_|«-^ (3.13) 

Since the transfer efficiency (3.9) depends strongly on the momentum transfer, this difference 
shows up in the ratio W+ /W- between the weights of the sidebands: 



— «exp(^.-J=exp(^.-3- 



exp ( TT— ) = exp ( TT — 3 ) (3-14) 



3.3. Quantum limit. — In the limit of small momenta Pi, Pf <C 1, one expects the evanescent 
wave mirror to produce the same resuts as a "hard" ideal mirror (Sect. 2.), since the wavelength 
of the incident atom is then longer than the characteristic decay length 1/2k of the mirror 
potential. In this limit expression (3.8) becomes: 

Wfi « ^e^PiPf (3.15) 

This result can be interpreted by noting that the vibrating evanescent wave mirror, in this 
regime [17], behaves as an ideal mirror moving as 

z^{t) = — ln(l + e smut) (3.16) 

(^) The asymmetry (3.13) is due to the non-linear dispersion relation Q. = hk^ /2M of de Broglie- waves. 
For light waves propagating in vacuum in one dimension, Q = ck, and this asymmetry vanishes. 
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For a small modulation amplitude e, this given a sinusoidal variation having amplitude zq = 

e/2K. Following the discussion of section 2, the weight of the first sideband is then given by 
I Ji(it)p !v u^/A where the argument u = ePi (using Eq. (2.4) and k = p/K). The flux into the 
sideband Pf is then given by (3.15), per unit flux in the incident state. 



4. Semiclassical path integral approach. 

In this section, we present a semi-classical perturbation method which allows us to derive the 
atomic wave function after reflection off the modulated mirror. This method is similar to the 
one developed for the problem of atomic wave diffraction by a thin phase grating [19]. It 
leads to a perturbed wave function which is phase-shifted with respect to the wave hmction 
obtained for a non-modulated mirror. The phase shift is simply the integral of the modulated 
potential along the classical unperturbed trajectory of the bouncing atom. In [19], this method 
was derived using a Fcynman path integral approach which allowed for a detailed study of 
the approximations involved in the derivation. For the sake of simplicity, we present here an 
alternative derivation based on a modification of the standard WKB treatment so that it can 
be applied to a time-modulated potential. 

4.1. Phase shift of the semiclassical wave function. — Consinder the wave function 
^^{z) which is an eigenstate of Hq at the energy E^. We look for a perturbed wave function 
of the following form 

^{z,t) = ^°{z)e-'^''l^ e^^{iSi{z,t)/n) (4.1) 

where Si (z, t) is a small correction introduced by the modulated part of the potential. The 
time-dependent Schrodinger equation gives: 

dSi f ih in d-'Si 1 (dSiV _ 

dt \m-^^{z) dz J dz 2M 9z2 2M \ dz J ' ^ ' ' 

In the semiclassical limit, where the de Broglie wavelenght of the particle is much smaller than 
the typical length scale of the potential Vo{z), we can approximate '^^{z) in the classically 
allowed region by the well-known WKB result: 

^^WKBi^) = exp (±1 f k{z')dz'^ (4.3) 

where the local wave vector k{z) is given by 

k{z) = ^^2M{Ei-Vo{z)) (4.4) 



Using now this result for '^'^{z), we can write 

- WKB 



±— ^ = v{z) (4.5) 



M^WKB^z) dz ' M 

where v{z) is the classical velocity of the particle moving in the unperturbed potential. It can 
be cither negative (before reflection) or positive (after reflection). Note that we have neglected 
in (4.5) the derivative of the prefactor of (4.3) which would lead to small corrections. The 
equation of evolution for Si can now be written 

dSi dSi in d^Si i /a^A' 

^ + = + 2M^ - 2M U7 



1884 



JOURNAL DE PHYSIQUE II 



N°ll 



This equation can be solved formally by the method of characteristics: using the characteristic 
curve Zc{t) defined by 

|.eW = K^c(i)) (4.7) 

the left hand side of (4.6) can be written 



d_ 
di 



v{z)i-]S, = ^S,{z,{t),t) 



dz 



dt 



(4.6 



The characteristic curve Zc{t) is the classical trajectory corresponding to a reflection in the 
unperturbed potential [17]: 



cosh^(t - to)/T 



(4.9) 



At the "bouncing time" t = to, the atom reaches its classical turning point. 
Substituting (4.8) in (4.6), we have the implicit solution 



Si{z,t) = - I dt'V,{z' = z,{t'),t') 

J — OO 

In the limit t —oo, (4.10) vanishes and (4.1) then reduces to the unperturbed wave function. 
We are interested here in the case t — to^^T, when the final time t is in the asymptotic region 
after the reflection. 

We now neglect the second term of the right hand side of (4.10). We will investigate the 
validity of this approximation in detail later on, but we can justify it here in a few words. We 
are neglecting a second derivative d^Si/dz and a second order term (dSi/dz)'^. The second 
derivative should have a small contribution in the semiclassical regime of interest here; its 
effect is mostly to correct the classical motion and to change the prefactor entering in (4.3). 
The second order term [dSi/dz)"^ should be small compared to the first order term entering 
in (4.2), since we expect Si itself to be a small correction to the unperturbed wave function. 

Using the expression (3.2) for Vi{z',t'), the integration in (4.10) now yields the phase shift 
[20]: 

Si{z,t) ^ _e /■* smojt' 
n h2M cosh2(i' - to)/T 

= -ePil3{(JT)sm(jto (4.11) 

Since we have t — to ^ t, this result does not explicitly depend on the upper bound t, up to 
small terms of order exp[— 2(t — to)/T]. 

We recover the function P{ujt) — (3{Q) defined in (3.10). The "bouncing time" to{z, t) of the 
classical trajectory (4.9) is fixed such that latter ends at time t at the position z: Zc{t) = z. 
Expanding the trajectory (4.9) for t — to ^ t, we find: 



(4.12) 
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where the "effective mirror" position ^eff equals [17] 



Vo_ 

4Ei 



We obtain finally 



Si{z,t) 

n 



I sin (u)t- |(2;-^eflF)) 



where the modulation index u is given by: 

u = ePiP{Q). 



(4.13) 
(4.14) 

(4.15) 



4.2. The final energy spectrum. — The reflected wave function in the asymptotic region 
KZ ^ 1 can now be written 



*fin(^,i) = C exp -{-Eit + piZ + hr] + Si{z,t)) 



(4.16) 



In this expression C is a normalization factor and // = r]{pi) is the phase shift of the wave 
function due to the reflection off the non- modulated evanescent potential. Replacing Si{z,t) 
by its expression (4.14), and using (2.5) to expand the result in terms of energy sidebands, we 
have 

^^niz,t) = Ce'" ane^^{-Ent+p„z) (4.17) 



where the n'th sideband has energy and momentum Pn{q = hojM/pi): 



En = Ei + nnw (4.18) 

Pn=Pi+nq (4-19) 

and its amplitude a„ equals: 

an = Jn{u)ex.p{-inq^es/^) (4.20) 

Note also that, if the phase of the modulated potential at t = is shifted by (j), then the phases 
of the sidebands are shifted accordingly: 

ujt^ujt-(l) an^ane'""^ (4.21) 

As usual in phase modulation problems, there are two limiting regimes for the result (4.17 
4.20), depending on the order of magnitude of the modulation index u. For a low modulation 
index, u ^ 1, the Bessel functions J„(u) have magnitude (u/2)l"l, and the diffracted spectrum 
consists essentially of the carrier n = and the two first sidebands n = ±1. In this regime, 
we recover the result (3.9) of the Born approximation for the first sideband, evaluated for a 
semiclassical momentum pi » Hk: 

7/2 1 

lair = \Mu)\' ^ ^ = -/P?P\Q) = (4.22) 

where we have assumed in addition Q <C Pi so that Q ~ AP. Note, however, that the 
semiclassical approach does not account for the asymmetry of the sideband weights (3.14); this 
property is related to the approximations involved (see below). 
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The semiclassical result (4.17) and (4.20) extends the Born result (3.9) into the region of 
high modulation index ?i ^ 1, where several lines arc present with an appreciable weight. The 
most intense lines correspond to n ~ ±u, and they lead to a velocity change for the atom 
Auat.max — uq/m. For e <C 1, where the vibration of -2m (i) is harmonic {2KZni{t) = e sin ivt), 
we can relate At;at,max to the "maximal velocity of the vibrating mirror" Umir,max = ew/2K: 

-max — ^'Ujnir.max 

(3{Q) (4.23) 

For a "hard" mirror formed with a potential step, one would expect simply Ai>at,max = 
2wmir.max- The reduction factor (3{Q) {P{Q) ^ 1 for Q > 1), which appears either in the 
modulation index u (see (4.15)) or in Ai;at,max (see (4.23)) is due to the "softness" of the 
potential. Indeed, the efficiency of the momentum transfer is proportional to the Fourier 
transform of the potential "seen" by the atom during the reflection (4.11): 

Voizcit)) = (4.24) 
cosh [t — toj/T 

the transform being evaluated at the modulation frequency w. Since this Fourier transform has 
a natural cutoff frequency of the order of t~^, the transfer efficiency decreases for modulation 
frequencies larger than this limit, i.e. Q > 1. The same reduction factor >6{Q) appears in the 
classical problem of a particle bouncing on a modulated exponential potential. The equation 
of motion for the particle is: 

= 2Kyoe-^''^(l + esinwt) (4.25) 
The energy change in the reflection can be written: 

dz d'^z 



f°° dz d^z 

Ei-Ei= M——^dt = -€Vouj cos{Ljt)e~^''^dt 
J —oo dt dt J —oo 



(4.26) 



using an integration by parts. We evaluate this last integral along the atomic trajectory in the 
non-modulated potential, obtaining 

El^El = -2eQ(3{Q) coswto (4.27) 
The maximal atomic velocity change (for coswto = ±1) is then: 

At;at,i 



Vi 

or 



= eO/3(g) (4.28) 



Ai;at,ma. = c-piQ) (4.29) 

K 

which is identical to (4.23). Note finally that exactly the same dependence on Q appears in 
the diffraction of atoms by a standing wave at oblique incidence [21, 19]. 

4.3. Validity of the semiclassical approach. — The first vahdity condition for a semi- 
classical approach requires 



Pi > 1 (4.30) 
i.e. an incident de Broglie wavelenght much smaller than the decay lenght of the potential. 
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A second constraint on the validity of (4.17) results from the use of the method of character- 
istics. Since the expression (4.11) for Si{z, t) has been obtained by integrating over the classical 
trajectory of the atoms in the absence of the modulated potential, we require that this classical 
trajectory is only slightly perturbed by the modulation. Otherwise the perturbative expansion 
underlying (4.11) would not be possible. This condition is fulfilled in two cases. One can first 
take a very small modulated potential (e <C 1); the frequency of the modulation can then be 
chosen freely. The other option is to take an arbitrarily large modulation factor (up to e = 1), 
but to impose a modulation frequency uj much greater than the characteristic "frequency" of 
the bouncing process l/r, where r is the reflection time. The modulated potential then induces 
a fast atomic micromotion, which is superimposed on the slow unperturbed bouncing motion. 
In the following, we focus on this second option, since it may lead to important transfer of 
momentum for the atoms. This condition can be written; 



WT = Q > 1 (4.31) 

In the previous calculation, an additional approximation is involved, which consists in ne- 
glecting in the expression (4.10) for Si{z,t) the contributions of the terms d'^Si/dz^ and 
{dSi/dzY- Before going further in a quantitative estimation of the corresponding error, we can 
point out two consequences of this approximation, which appear clearly in the result (4.17). 
First, (4.17) does not strictly fulfill the dispersion relation for matter waves = p^/2M. 
Equations (4.18) and (4.19) only constitute a linearized version of this dispersion relation for 
small nq, and this induces a small phase error, that should be compensated by the contribu- 
tion of these two partial derivatives neglected in (4.10). Secondly the momentum distribution 
deduced from (4.17) is a symmetric comb centered on pi, with sidebands whose weight is pro- 
portional to |a„p = |a_np. However we know from the Born treatment (Sect. 3) that an 
asymmetry appears in the spectrum, when the parameter / Pi becomes of the order of 1 or 
larger (see (3.13)). We show now how one can recover this validity condition from a detailed 
analysis of the contributions of the two partial derivatives mentioned above. 

We estimate the magnitude of the second term of the right hand side of (4.10), in order to 
determine the region of parameter space where it can be neglected with respect to the leading 
term. As an estimation we take Si{z,t) = before the boimcing time (t < to)j and we use for 
t > to the asymptotic expression for s[^\z,t) given in (4.14) (^). We obtain in this way, for 
t > to: 



dS 

—— = uqcos{u)tQ) (4.32) 
oz 

|^ = ^sin(c.t„) (4.33) 

We now choose a time t such that t — to is larger than the reflection time t and we evaluate 
the two contributions that have been neglected in the integral appearing in (4.10): 



1 

2Mh 



2M 

t 

-oc 



dz'^ 



{z',t')dt' 



Q'^ t-tp 

, t ~ to 
"P, 2t 



(4.34) 
(4.35) 



(^) A more precise evaluation of Si would require lengthy calculations on this semi-classical framework, 
due to spurious divergences of the WKB method around the classical turning point. 
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The relative magnitude of these two terms depends on the index of modulation. Equation 

(4.34) is the leading term in the regime of low modulation index {u <C 1), whereas (4.35) is 
dominant in the high modulation domain (m 3> 1). 

We note that the magnitude of the two terms (4.34) and (4.35) increases linearly with time, as 
expected for corrections to the phase of ^'fin(z,t) whose role is to restore the correct dipersion 
relation for the de Broglie wave. In practice, we can evaluate these two terms for a time t 
located well after the classical turning time to. We take for instance: 

t - io - 4r (4.36) 

After the time t, the values of the a„ coefficients appearing in the final wave function (4.17) 
do not change anymore, and we can impose "by hand" the correct dispersion relation; in other 
words, we then replace the wave function (4.17) by 

^U^, t) = Ce'^ J + P'n^) (4-37) 

n 

with = y/2MEn. 




Fig. 2. — Regions of validity of the semiclassical approach, as a function of Q and Pi = pi/hn. The 
modulation depth was taken as e = 1. The semiclassical approach is valid in the region between 
the lines (b) and (d). The line (c) [u — 1) separates the regimes of high (above) and low (below) 
modulation index u. (a) Line (3 = 3, (Eq. (4.30)) defining the region for a fast classical micromotion 
of the atom, (b) Limit 2{APmax)''^ = Pi (Eq. (4.39)). (d) Limit 2Q^ = Pi (Eq. (4.38)). (c) Limiting 
condition pi = Hk for the semiclassical regime. Crosses: parameters corresponding to the numerical 
solution (Fig. 1 and Sect. 5). 



In the regime of low modulation index, the leading correction (4.34) is small compared with 
the main contribution to Si/h if: 

2Q2 < p. (4.38) 

As expected, we recover here the condition required for having a quasi-symmetric spectrum. 

In the regime of high modulation index, the main contribution to the phase factor Si/h is 
very large compared to 1, and we now have to require that the leading phase correction (4.35) 
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be small compared to 1. We then obtain: 



2(AP^ax)' « Pi 



(4.39) 



where APmax represent the maximal appreciable momentum transfer, measured in units of Hk: 



The validity region delimited by the two equations (4.38) and (4.39), together with (4.30) and 
(4.31), is plotted in figure 2. 



5. Reflection of wave packets. 

In order to check the approximate results (3.8) and (4.20), we have performed a direct numerical 
integration of the Schrodinger equation, giving the evolution of an atomic wave packet bouncing 
on the atomic mirror. As in sections 3 and 4, we restricted ourselves to the one-dimensional 
problem, since the atomic degrees of freedom parallel to the plane of the mirror factorizc out. 
We choose an initial Gaussian wave packet and we integrate the evolution of the atomic wave 
function in position space using a 4'th order Runge-Kutta algorithm [22]. 

The position of the initial wave packet is chosen far enough from the mirror position so that 
its propagation towards the mirror is initially the same as for a free particle. For instance, 
for the evolution shown in figure 3, the initial wave packet was centered at KZi = 13, with a 
standard deviation nSzi = 2. The initial momentum was pi = 100?i«:, with a standard deviation 
Spi = 0.25 Hk. We have checked that the final momentum distribution is independent of the 
value of Spi, provided that Spi <^ q, i.e. in the limit of well-resolved orders. 

In our calculation, the particles are confined in a square box whose limits -25 
and KZmax = 25. The boundary conditions for the wave function V'(-z) are ■(/'(•^min) = i^{zma,x) = 
0. The expression of the potential in the ^ > domain is the same as the one used in section 
3 (sec Eq. (3.1) and (3.2)). In the z < domain, the potential is chosen to be null. This 
form for the potential allows us to study the fraction of atoms that reach the point z = 0. 
In a real experiment these atoms would actually hit the surface of the dielectric supporting 
the evanescent light wave. They would either stick to the dielectric, or be rccmittcd with a 
thermal velocity; in any case they would be lost for the subsequent use of the bouncing atomic 
beam. For the values of the parameters used in the following examples, we have checked that 
this fraction of atoms is always negligible, even for a 100% modulated potential (e = 1) (''). 

After the reflection, i.e. when the reflected wave is in a region where the influence of the 
potential is negligible, we calculate the momentum distribution \-ijj{p)\^, where 'ip{p) is the 
Fourier transform of tpiz)- Two example momentum distributions are shown in figure 4. We 
have checked that the result for |?/j(p)p does not depend on the time at which the Fourier 
transform is taken, as expected from the simple evalotion of ■)/'(}?) after reflection: 'ip{p,t2) = 
'0(p, <i) exp(zp^(t2 — ti)/2Mfi). On the other hand, the shape of the position distribution 
|V'(2)p changes long after the reflection. For instance the oscillations appearing in figure 3, 
which result from the superposition of wave packets with momenta pi, pi ± q, will eventually 
disappear when the three wave packets become separated from each other because of their 
different group velocities. 



AP^ax =uQ = eP,Qf3{Q) 



(4.40) 



(^) For the practical design of an experiment, one should include, for safety, the van der Waals 
attractive potential in order to estimate more precisely the fraction of sticking atoms. 
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momenturn (p/ft/c) b) 

Fig. 3. — a) Reflection of a wave packet at the modulated mirror. The position distribution at t = 
has been cut for the sake of clarity, its actual maximum height is about 15 times the height of the 
incident wave packet. The parameters used arc: initial position nz^ = 13 with standard deviation 
Skzi = 2, initial momentum p\ = 100 Tin with standard deviation 5pi = 0.25 Tin, modulation frequency 
w = 5/r so that 9 = 5 Tin, modulation depth e = 1. b) Final momentum distribution. 



We now compare the predictions of this numerical treatment with those of the approach 
presented in section 4. To this purpose, we consider again the initial Gaussian wave packet: 

^(p,0) = exp (-^^^) exp . (5.1) 
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Fig. 4. — Final momentum distributions, given by the numerical solution (solid line) and the semi- 
classical approach (dashed line). The parameters used are: initial momentum p\ = 1007i/t, momentum 
width 5pi = 0.25 Uk (standard deviation), modulation frequency ui = 4.2/r (momentum transfer 
Ap = 4.2 hn). The situation adopted in (a) (modulation depth e — 0.6) maximizes the product of the 
carrier and first sideband; that in (b)(e = 1) ma^ximizes the product of the sidebands n = ±1. 



Using the approximate results of section 4, the final wave packet is given by 

= E-e-exp (-m^) «P - li^i) (5., 

where the a„ coefficients, given in (4.20), must be evaluated a.t p — nq. The results of the 
two approaches are shown in figure 4, the solid line for the numerical solution and the dotted 
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line for the semiclassical approximation. The initial wave packet is the same as above. The 

momentum transfer is q = 4.2 Hk and the modulation amplitude e equals 0.6 (Fig. 4a) and 
1 (Fig. 4b). We see that the agreement between the predictions of the two methods is good, 
although not perfect. In particular, the numerical treatment shows an asymmetry between 
the heights of the two sidebands \ijj{pi ± q)\'^ , while the approximate treatment predicts equal 
sideband weights, because of the relation | Jra(a;)| = | J_„(a;)| (see discussion in Sect. 4). 

A more systematic comparison of the predictions of this numerical approach with the results 
derived with the Born approximation and with the semiclassical approach is presented in figure 
1. We have determined, for the same initial wave packet as above, and for several values for q, 
the heights of the sidebands n = ±1. We see that the agreement between the three methods, 
in their expected range of validity, is quite satisfactory. 

6. An atom interferometer. 

The vibrating mirror can be used in a number of ways to make an interferometer. In what 
follows we will consider the case of atoms normally incident on the mirror surface, and making 
three or more bounces on it. In other words, wc rely on the possibility of having a source of 
slow atoms released above the mirror [23, 24]. With a fast beam of atom, either which could 
not be reflected at normal incidence, or for which the time to perform repeated bounces is too 
long, one would need several mirrors — used at grazing incidence if necessary. 

Figure 5 illustrates an interferometer based on three consecutive bounces on a vibrating 
mirror. This is similar in conception to an interferometer constructed from three diffraction 
gratings [6, 25]. For an incident energy Ei, we consider two output channels with energies 
equal to and E^ + hu). Each of these two channels can be reached via two paths, which are 
shown by full lines in flgure 5; the other paths, shown dashed in figure 5, do not contribute 
since we assume that the mirror is "turned off" when these paths hit it. 




A B D C 

r\AAA/^ rww^i rww\_ 

sin{ut - tp,) sLn(wf - (p^) sin(u t - ^p,) 



Fig. 5. — Space-time diagram showing a monochromatic interferometer using three bounces on a 
vibrating mirror. Tiio paths siiown lie at local minima of tiic classical action. At each reflection, paths 
corresponding to sidebands of order higher than 1 are omitted for clarity, and paths which return to 
the mirror when the latter is non-reflecting are shown dashed. 

The probability amplitude to exit in a given channel is calculated using the propagators for 
the region where the atom is in free fall, and the phase shifts corresponding to the interaction 
with the mirror. The free fall propagator can be evaluated using the integral of the action 
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along the relevant classical paths. We write these in the form exp(ia^B) for the various paths 
as shown in figure 5. Wc write (f>m the phase of the mirror vibration at the m**^ bounce. 
This phase is dictated purely by the origin in time of the sinusoidal modulation of the mirror 
for that bounce. Using the expressions (4.20) and (4.21) for the coefficients a„ giving the 
amplitude for transmission into the n'th sideband, we now get the two probability amplitudes 
for transmission into the channels: 



A{Ei + Two) = 



+ al^^aL14^^e*('^i-'*^)c'("^"+«°'^) (6.1) 



+ a^^^aL^^ap^e'('^i-'^=^+'^^)e'("-*°+"°^) (6.2) 



where the superscript (m) indicates the bounce number. For simplicity we have omitted here 
the contribution of the phases ri{p) which appear in (4.17). This contribution is the same for 
the two paths and it cancels out in the final interference pattern. 

By choosing a symmetric geometry, we impose the conditions aAB = cxdc and a ad = cxbc 
so that the final interferometric phase is given simply by the refiections; it can be written 
9 = (pi — 2(j)2 + (j)3- Note that the other phases associated with the reflections, g^cff and, 
as already noted, r?(p), exactly cancel out because of the symmetry. In other words all the 
phases which depend on the incident atomic momentum disappear, and the interferometer 
will produce high-contrast fringes even when illuminated by a beam of atoms having a broad 
momentum distribution. The "tradional" three-grating interferometer has the same property. 
Note also that, for any phase grating including the present one, the a„ coeflacients satisfy 
aQtt-i = — aoa|. This ensures that the total probability for ending in one channel or the other, 
[^(i^i)!^ + |>l(-Bi + 7icj)p, is independent of the interferometric phase 6. 

The fringe amplitude for this interferometer can be written 



^= 2Re 



(aoat)(i)(aia*_i)(2)(a_iaS)(3)e'*' (6.3) 



This quantity would give = cos 6 for an ideal interferometer such as a Mach Zehnder. Using 
the semiclassical approximation (4.20), and denoting Um the index of modulation for the m'th 
bounce, one has 

T = 2iMui)Mui)){Ji{u2)fiMu3)Jiiu3))cos0 (6.4) 

This fringe amplitude JF is optimised for ui = U3 = 1.08 and U2 = 1.84, and we get |J-"|max = 
0.078. To maximize in these conditions the separation of the interferometer arms we choose 
e2 = 1 at the middle boimcc, and (4.15) then gives ci = £3 ~ 1.08/1.84 ~ 0.6 for the first 
and third bounces. For example, if the initial momentum is pi = 100 Hk, (4.15) leads to 
Ap = 4.2 Hk for these values of the e's and u's. This gives a good indication of the optimum 
case: the numerical solution of the Schrodinger equation gives the maximum |.?'|max — 0.081 
at Ap ~ 4.2 hn, for the same value of pi (see Figs. 4a and 4b). To optimise instead the fringe 
contrast, one would use the case |oo| = |ai| for the first and third bounces, producing 100% 
contrast, with a slightly reduced fringe amplitude. 

As an example, consider a cesium atom normally incident on an evanescent wave near- 
resonant with the atomic transition 6S1/2 — * QP3/2- Taking k ~ 27r/852 nm, the momentum 
100 hn corresponds to a velocity of 35 cm/s, and a boimcc height of 6 mm, which is readily 
realisable in practice. The momentum change 4.2 ?ik is produced with w = 27r x 1.7 MHz. 
The two interferometer arms are separated in distance by about 0.5 mm, when both are near 
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Fig. 6. — Interferometer design in which the interference phase is sensitive to the acceleration due to 
gravity. Since the vibrating mirror can only transfer a few Uk of momentum, a practical design would 
involve several bounces before the arms are recombined, as described in the text. 



their maximum height, and in time by about 3 ms (the difference in arrival times at the second 
bounce). Thus macroscopic path separation can be obtained. 

The interference phase will be destroyed not only by mechanical fluctuations of the dielectric 
surface supporting the evanescent wave, but also by intensity fluctuations of the evanescent 
wave itself. To estimate the effect of the latter, we note that if the intensity fluctuates by a 
fraction e, then the effective mirror position (Eq. (4.13)) moves by e/2K. The change produced 
in the interferometer phase is then tp/fin radians; this can be understood simply as a higher 
sensitivity to fluctuations when the de Broglie wavelenght is small, or, equivalently, when the 
interferometer area is large. One must eliminate not only fluctuations in time but also in 
the transverse spatial profile of the laser beam forming the evanescent wave "mirror" . To 
make a flat mirror, one could use a Gaussian laser beam of waist w internally reflected from a 
concave glass surface of radius of curvature R: when w = ^J2R/ k the glass surface curvature 
compensates the Gaussian fall-off of the laser beam intensity profile, and the mirror presents 
a flat reflecting surface for the atoms. A transverse magnetic quadrupole field could help to 
confine the atom without perturbing the vertical motions of the interferometer arms. 

The most obvious first use of the vibrating mirror interferometer is simply to investigate the 
preservation of coherence during a reflection off the evanescent wave mirror — this could be 
done by introducing more and more bounces between the splitting and recombining of the two 
interferometer arms. A second use is as a probe of the acceleration due to gravity. For this, 
one would not use the symmetric arrangement which we have considered so far, but instead 
one allows the two arms to bounce many times until they arc recombined "automatically" after 
m (and m + 1) bounces, where mp' = (m + l)p, with p'^ — p^ + 2Mhuj (Fig. 6). The mirror 
is vibrated only for the initial and final bounces; in between it is stationary. For this case, the 
interferometer phase is dominated by the contributions from the free flights between bounces. 
One finds that while such an interferometer provides a sensitive probe of the acceleration due 
to gravity, it is also highly chromatic, producing fringes only for a very small class of incident 
momenta centred around p. 

Conclusion. 

To summarize, we have presented here the principles of a vibrating mirror for atoms, which 
constitutes a novel beam splitter for atomic de Broglie waves. We have investigated two 
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analjdiical theoretical approaches for this new scheme, in order to determine the amphtudes 
and the phases of the outgoing atomic waves. The results are in excellent agreement with a 
numerical approach to the problem. 

The typical momentum transferred by a vibrating mirror formed with an evanescent light 
wave is a few photon momenta. Although this is not as large as the transfers obtained by some 
other devices, we believe that it should provide, because of its conceptual simplicity, a conve- 
nient tool for atom optics and interferometry. Indeed, unlike other atom optics components 
such as micro-fabricated gratings, it is quite easy to change rapidly (on the microsecond scale) 
the modulation factor of the light wave forming the mirror. This gives in return a direct con- 
trol upon the phases and intensities of the diffracted de Broglie waves, and this allows one to 
conceive simple and useful atomic interferometric devices, such as the ones shown in figures 5 
and 6. 
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